Dynamic Response Functions from Angle Resolved Photoemission Spectra 
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We introduce a formalism for calculating dynamic response functions using experimental single- 
particle Green's functions. As an illustration of this procedure we estimate the dynamic spin response 
of the cuprate superconductor Bi2Sr2CaCu2 0s+i. We find good agreement with superconducting 
state neutron data, in particular the (ji, n) resonance with its unusual 'hourglass' shaped dispersion. 
We anticipate our formalism will also be useful in interpreting results from other spectroscopies, 
such as optical and Raman responses. 

PACS numbers: 74.25.Jb, 74.25.Ha, 74.72.Hs, 79.60.Bm 
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The linear response to an external probe as a func- 
tion of momentum and frequency is of great importance 
in elucidating the properties of complex materials. Ex- 
amples include various two-particle correlation functions 
involving spin, current and charge as measured by inelas- 
tic neutron scattering (INS), nuclear magnetic resonance 
(NMR), optical conductivity, and Raman scattering ex- 
periments. On the other hand, angle resolved photoemis- 
sion spectroscopy (ARPES) [l[ directly gives information 
about single-particle excitations of a system. The re- 
sponse function of a system can be expressed in terms 
of a two-particle correlation function of the observable 
to which the external probe couples. The goal of this 
paper is to develop an approach to use single-particle 
spectroscopy data to gain insight into two-particle corre- 
lation functions. In particular we focus here on using the 
Green's functions obtained from superconducting state 
ARPES data in the high T c cuprates to compute the dy- 
namic spin susceptibility, which we then compare with 
INS data 0. 

From a theoretical point of view, dynamic response 
functions are difficult to calculate in general and many 
different approximate formalisms exist in the literature. 
For instance, there are two rather different approaches 
for computing the dynamic spin response for the high 
T c cuprate superconductors. The first is based on the 
random phase approximation (RPA) [3| and related di- 
agrammatic formulations [J]. This approach not only 
assumes momentum is a good quantum number, but also 
that the spin and charge degrees of freedom are cou- 
pled. The second is based on spin ladders separated by 
one-dimensional domain walls known as stripes. In this 
formalism spatial inhomogeneity is important, and the 



charge sector is assumed to be secondary when calcu- 
lating the spin response [H]. Despite the quite different 
physics underlying these two schemes, the results for the 
calculated spin response function of the cuprates are sim- 
ilar - one of the current dilemmas facing the field of high 
T c superconductivity. It is thus important to go beyond 
a purely theoretical approach and directly employ infor- 
mation obtained from one experiment (ARPES) to make 
progress on interpreting the dynamic susceptibility mea- 
sured by another (INS). 

We use a formalism based on a diagrammatic k-space 
approach which goes beyond RPA in that it uses fully 
dressed Green's function obtained from ARPES data on 
Bi 2 Sr2CaCu208+(5 (Bi2212). We compare the calculated 
superconducting state susceptibility with INS data. We 
obtain the (tt,tt) resonance seen in many cuprates Q, 
including Bi2212 Q, and also its unusual 'hourglass' 
shaped dispersion as observed in YBa2Cua07_5 (YBCO) 
@>[H and more recently in Bi2212 [9]. We also find that 
the magnetic dispersion is sensitive to the momentum 
dependence of the effective interaction used to calculate 
the susceptibility 

We use ARPES spectra from a near-optimal sample 
(T r ~ 90K) of Bi2212, the data having been presented pre- 
viously [la [ll|. While a resonance peak was observed 
in this material some time ago [6] , a more detailed study 
with results similar to the much more extensive INS data 
for YBCO, has appeared only recently 

Quite generally two-particle correlation functions can 
be written in terms of single-particle Green's functions 
and vertex parts [l2| . The lowest order term contributing 
to the spin susceptibility (the bare polarization bubble) 
in the superconducting state can be written as [l3| 
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where Im denotes the imaginary part of the normal and 
anomalous Green's functions G and F , and Xo an d Xo 
denote the GG and FF contributions to xo respectively. 

We next describe in detail how ImG is extracted 
from ARPES data and return later to the question of 
estimating the contribution of InxF (which is not di- 
rectly measured). ARPES probes the occupied part of 
the spectral function leading to the intensity Z(k, u>) oc 
np{u>) ImG(k, ui), where np{uj) is the Fermi function [14j |. 
In order to extract ImG from raw data we need to address 
several issues including data normalization, background 
subtraction, and removing the effects of the Fermi func- 
tion. In addition we need to extend ImG to ui > to 
calculate xo- 

Starting from raw data, we first subtract the constant 
signal at uj > (due to second order light). Next an 'un- 
occupied' state spectrum at a k far from kp is used as an 
energy-dependent background [l5j |. The subtraction is 
performed by normalizing the background to each spec- 
tra at a given binding energy, ui c (320 meV for the data 
set in question), and then subtracting it [lj|. This ef- 
fective spectral function represents the the renormalized 
band near the Fermi energy. Finally, we divide the data 
by a resolution 17J broadened Fermi function to obtain 
ImG(k,w) for ui < 0. 

The next step is to determine the unoccupied part of 
the spectral function, ImG(k, ui) for ui > 0, which cannot 
be obtained directly from ARPES data. We obtain this 
by invoking particle-hole symmetry with respect to the 
Fermi surface: k^, ImG(k^ +k, ui) — ImG(k^ — k, — ui), 
where k is directed along the normal to the Fermi surface. 
This assumption should be reasonable in the supercon- 
ducting state of optimally doped cuprates over an energy 
range in excess of the gap, as evidenced by the approx- 
imate particle-hole symmetry seen in tunneling experi- 
ments |l8j. We have also checked that this assumption 
does not qualitatively affect our final results for xo (by 
using ImG with p-h asymmetry put in by hand). Wc 
then normalize the obtained ImG so that the integral of 
the spectral function (— ImG/7r) is equal to unity over the 
energy range of ±ui c . This minimizes the effect of dipole 
matrix elements. Now we may use the ImG derived from 
ARPES to calculate Xo ( we w iU discuss Xo later). Fi- 
nally, to perform the k-sum in Eq. ([1]), the ARPES data 
are interpolated to a regular grid and then reflected using 
square lattice group operations to fill the first Brillouin 
zone [llj]. We used a 100 x 100 grid. 

Fig. la shows the calculated Imxff at T=40K (super- 
conducting state) as a function of the momentum transfer 




40 80 120 160 

a (meV) 



40 80 120 160 

a. (meV) 



FIG. 1: (a) Imx§ and (b) Rexo in the superconducting state 
as a function of frequency for momenta along q = 7](tt,tt). 
The various curves, labeled by r\, are offset for clarity in panel 
(a), noting that Imxo = at = 0. 



q along the zone diagonal. We note that Imxcf is greatly 
suppressed at low energies due to the gap to particle-hole 
(p-h) excitations in the superconducting state and then 
increases quite sharply. The p-h gap is in general given 
by the sum of the superconducting gaps at two points on 
the Fermi surface separated by the wavevector q. The 
Q = (it, it) vector connects the hot spots (ek = ek+Q=0) 
which are not too far from the zone boundary in Bi2212, 
and thus the hot-spot gap is comparable to the one at 
the antinode. Consequentially at Q = (it, it), we see in 
Fig. la a large gap whose midpoint is around 80 meV, 
roughly twice the maximum d-wave superconducting gap 
of ~ 40 meV at the antinode. We note that the thresh- 
old is quite broad (~ 40 meV) as a result of the intrinsic 
broadening of ImG arising from self-energy effects as well 
as resolution. As q decreases from Q, one sees the p-h 
gap decrease due to the d-wave anisotropy of the gap, 
and then disappear at q n ~ (0.76, 0.76)7r. q„ is the 
wavevector corresponding to node-node scattering with 
the d-wave gap vanishing at the nodes. For q < q n , the 
p-h gap reappears 0]. 

In Fig. lb we show Rex^ obtained from Eq. (TTJ) [lj|. 
First concentrating on Q, we note the presence of a 
peak that corresponds to the gap midpoint of Fig. la 
as expected from Kramers-Kronig relations. This peak 
is broadest for q n where the p-h gap vanishes in the imag- 
inary part. 

We now turn to Xo ■ Since ImF is not available from 
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experiment, we estimate the FF term as follows. We 
calculate the BCS Xcxbcs and Xo^bcs irom Eq. |T|) us- 
ing the bare BCS Green's functions Go(k, uo) = (uj + 
£k )/(c 2 -el- A 2 ) and F (k,u) = A k /(^ 2 - j - A 2 .) 
with the experimentally measured dispersion [20| ek and 
the measured Ak, which we find to be proportional to 
(cosk x — cos k y ) for this sample. We define the ratio 
of the real and imaginary parts, given by a^(q, 51) — 
R" e XaBCs/R- e XoBCS anc ^ similarly f° r ^(^ft)- We then 
assume that the missing contribution Xo ma y be ac- 
counted for by Rexo = R e Xcf + ReXo" — (1 + an)Rcxo 
and Imxo = ImXo' + ^ m Xo — (1 + a/)fniX(?- We will 
discuss below the extent to which our final results are 
affected by this approximation. 

In order to carry out comparisons with INS data, or 
other probes such as NMR, the full spin susceptibility x 
is needed. The most common approximation is to use 
the RPA form [3] 



x(q,ft) 



xo(q, ft) 



l-tf(q)Xo(q,n) 



(2) 



where U is an effective interaction. In this paper, we will 
assume two limiting forms for this effective interaction, 
one where it is a constant (Uo), the other where it has 
the form U(q) = —Uo(cosq x +cosq y )/2 corresponding to 
superexchange between near neighbor copper sites. 

The experimental results for xo presented in Fig. I are 
similar to those obtained from BCS theory [3j especially 
at low energies where the incoherent spectral weight in 
the ARPES data is small. Within BCS theory, which 
uses bare Green's functions, one has a true gap in Imxo 
at Q and a corresponding log divergence in Rexo ■ These 
features still persist in Fig. I albeit broadened by self- 
energy (and resolution) effects. As such, for some fre- 
quency smaller than the gap, one will obtain a pole in 
X when I — [/(q)Rcxo(q, ft) = provided Imxo is small 
at the frequency of interest. This pole represents a col- 
lective mode, known as the spin resonance at Q, which 
is prominently observed in INS data for YBCO Q and 
Bi22I2 Following this logic, we fix U by fit- 

ting the energy (44 meV) of the INS spin resonance at Q 
for optimal doped Bi22I2 @]. 

In order to compare the results of our approach with 
INS data we plot the imaginary part of the full suscep- 
tibility x as calculated from ARPES data as discussed 
above for constant q cuts in Fig. [5] and constant f2 cuts 
in Fig. [3J The left panels assume a constant U, the right 
panels the near neighbor exchange form E/(q). Let us 
first consider Fig. O For constant U (left panel), the 
resonance traces out a pronounced downwards disper- 
sion which terminates at the node-node scattering vec- 
tor q ra , as seen in the INS data. For q < q n , a distinct 
second branch appears, which is broad and weak, that 
disperses upwards. The change of behavior at q„ cor- 
responds to the so-called silent band effect and second 
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FIG. 2: Imx as a function of energy for several momenta 
q = rj(n,n) assuming (a) a constant U — Uo or (b) a near 
neighbor exchange U = U(q). The curves, ranging from 
?7=0.60 (top) to 77=1 (bottom), are offset for clarity. Some 
curves are rescaled for visual clarity. 



mode mentioned in connection with INS data of over- 
doped YBCO [53, SH and Bi22I2 [§]. But interestingly, 
the upper branch of the 'hourglass' is not apparent. In 



contrast, for U(q) (right panel), the mode is almost dis- 
persionless near (7r,7r), then shows an upward dispersion 
for momenta q < 0.9(tt,tt). This difference is a direct 
consequence of the relatively weak momentum depen- 
dence of Rexo (Fig. lb) coupled to the decrease in U(q) 
away from Q — (tt, tt). 

The behavior of the dispersion observed in constant q 
scans can be contrasted with that from constant Q scans 
(Fig. [3|). For constant U, both types of scans yield a 
qualitatively similar mode dispersion (solid dots in Fig. 
[3ji). However, for U(q), the dispersion obtained from the 
constant fl cuts is hourglass-like, with an upward and 
downward branch that merge at (it, tt), in good agree- 
ment with INS data in underdoped YBCO [8]. We note 
that this downward branch is not visible in the constant 
q cuts in Fig. [2b. This difference is analogous to the 
different ARPES dispersions that one finds from energy 
distribution curves as compared to momentum distribu- 
tion curves. 

We have also generated results involving only the GG 
contribution by setting a = 0. In order to compare 
with the results having both FF and GG contribution, 
we rescale Uo to maintain the same resonance energy at 
Q. Only minor differences are found between this and 
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FIG. 3: Imx as a function of momentum q = rj(n,n) for 
several energies assuming (a) a constant U = Uo or (b) a 
near neighbor exchange U = £/(q). The curves, ranging from 
51=68 meV (top) to 20 meV (bottom), are offset for clarity. 
Some curves are rescaled for visual clarity. 



the full calculation which includes both GG and FF con- 
tributions. This lack of sensitivity to the inclusion of 
FF terms is a consequence of the d-wave symmetry of 
the superconducting order parameter where Xo an d Xo* 
reinforce one another near (tt,tt). (In contrast, for an 
s-wave superconductor, Xo an d Xo are opposite in sign, 
and there is no spin resonance.) 

Comparing our results with the INS data and ear- 
lier RPA calculations, we arrive at the following con- 
clusions. First, the self-energy effects present in the 
ARPES-derived Green's function do not affect the low 
energy physics of spin excitations, such as the existence 
and sharpness of the (it, tt) resonance, or the mode dis- 
persion. Second, vertex corrections do not play a major 
role in the spin channel, except possibly in the overall 
scale of U. Third, the magnetic dispersion is sensitive 
to the q dependence of U. Finally, our results provide 
strong evidence for the interpretation of the resonance 
peak as a spin exciton 0]. 

To summarize, we computed the polarization bubble 
Xo using experimental Green's functions derived from 
ARPES spectra, and the full dynamic spin response ob- 
tained from a diagrammatic formalism assuming either 
a constant or a near neighbor exchange interaction. Al- 
though this analysis requires several approximations, we 
find surprisingly good agreement with inelastic neutron 
scattering data for high temperature cuprate supercon- 
ductors. Our results demonstrate a close relation be- 



tween experiments probing the spin and single-particle 
excitations. Our formalism is quite general, and can be 
used as well to compute other response functions, such as 
the current-current response function measured by con- 
ductivity. 
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